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Abstract 

The stability is analyzed of the magnetic junction coUinear configura- 
tions against small fluctuations under amplitude-modulated current with 
CPP mode. High spin injection is assumed. Under parametric resonance 
conditions, with the modulation frequency twice the precession frequency, 
instability is possible of one, or another, or both the coUinear configura- 
tions. When the dc component of the current density exceeds the instabil- 
ity threshold of the antiparallel configuration, the parametric instability is 
suppressed by nonparametric one which is induced by the dc current. The 
parametric instability manifests itself as lowering the threshold of the dc 
current density in presence of the high-frequency current, such an effect 
has been observed in experiments repeatedly. 

Studying the magnetic junctions under spin-polarized current showed 
a number of interesting phenomena, such as the switching of magnetic con- 
figuration [T], spin- wave generation [2], etc. Such effects may appear in 
nanoscale level, because the corresponding characteristic lengths, namely, 
the exchange length and the spin diffusion length, are of order of 10 nm [3]. 
This allows using the effect for high-density information recording by elec- 
tric current, which is unattainable for switching magnetic elements by 
magnetic field alone. 

The current-driven switching magnetic junctions is accompanied often 
with the magnetization oscillation and the other high-frequency effects 
(see, e.g., [2] |31 El [S])- In this connection, an interesting problem occurs, 
namely, the influence of spin-polarized high-frequency current on magnetic 
junctions. 

The microwave effect on the switching magnetic structures by spin- 
polarized current was studied in Refs. [7ll8ll9l [T0llll| . Lowering threshold 
value of dc current density was found in presence of a high-frequency 
current. Nowadays, interpretation of this effect is restricted mainly with 
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Figure 1; Scheme of the magnetic junction. 



qualitative arguments and numerical simulations without developing any 
consistent theory. Without having pretensions to constructing a compre- 
hensive theory, we should like to pay attention to a possible mechanism 
of the effect related with occurring instability of the magnetic junction 
configuration under parametric resonance conditions. 

Let us consider a magnetic junction consisting of a pinned ferromag- 
netic layer 1, free ferromagnetic layer 2, and nonmagnetic layer 3 closing 
the electric circuit. Amplitude-modulated electric current with density 



flows perpendicular to the layers. 

We assume j < j, so that the total current flows always in the same 
direction, corresponding to the conduction electron flux in the 1 — ;> 2 3 
direction. We take, also, fir ^ 1, r being the electron spin relaxation 
time, that allows describing the alternative current effect by the same 
equations as the direct current one with changing constant current density 
by time-dependent one, j{t). 

We take a coordinate system with x axis perpendicular to the layers, 
yz plane parallel to layers, and x coordinate counting off the interface 
between 1 and 2 layers (Fig.[T|. The pinned layer magnetization vector Mi 
is directed along z axis. The free layer magnetization vector M originally 
(without current) is collinear (parallel or antiparallel) to the pinned layer 
magnetization vector Mi — {0, 0, il/i}. We are to investigate stability of 
such a state in presence of the current. 

The free layer thickness L is assumed to be small compared to the 
spin diffusion length in that layer and the domain wall thickness. In 



= j +jcosnt 



(1) 



2 



such a case, the macrospin approximation is valid, in which the free layer 
magnetization is taken as spatially uniform, while the presence of the 
current is taken into account by means of additional terms proportional 
to the current in the equations of the motion for the lattice magnetization 
(see Refs. [I1[T3] for detail). 

In the linear approximation on small deviations of the M vector from 
the coUinear position, the Landau-Lifshitz equations take the form 

dKLj, , — dMy 



+ f,M,—^ + 'ylH, + HaM 
at at \ 

+ P(t)My + K{t)TUMa: = 0, (2) 



-P{t)M^ + K{t)M,My = 0. (3) 

Here M = M/|Mj is the unit vector along the free layer magnetization, 

Mz — ±1 is the original stationary coUinear position of the vector, Mx,y 
are the small deviations from that position, H is the external magnetic 
field. Ha is the anisotropy field of the free layer, k is the Gilbert damping 
constant (k <C 1 is assumed), 7 is the gyromagnetic ratio. We assume 
\Hz I < Ha that excludes possibility of switching the junction by magnetic 
field without current. 

The functions K{t) and P{t) describe the spin-polarized current effect 
for two different interaction mechanisms between the conduction electrons 
and magnetic lattice, namely, spin transfer torque [14II15| and appearance 
of a region in the free layer with nonequilibrium electron spin polariza- 
tion [16M17| : for brevity, we name these mechanisms as torque and injec- 
tion ones, respectively. Note that the contribution of the injection mech- 
anism in the Eqs. ((2]) and ((3]) has the same form as the magnetic field, so 
that one may say about some effective magnetic field H^a = Hz + P{t)/^. 
Note, however, that unlike the similar effect of the parametric instabil- 
ity under longitudinal pumping by high-frequency magnetic field [18], the 
spin-polarized current effect, as mentioned above, is of local nature with 
characteristic scale of order of the spin diffusion length I 10 ~ 100 nm. 

Further we consider the case of high injection when the following in- 
equalities fulfill: 

Z2/Z1, Z2/Z3 <:\<^1, (4) 

where 

^^ = ^7(1^ (. = 1,2,3), (5) 

h is the spin diffusion length for ith layer, Qi = (a^^ — (Ji_i)/ai is the 
conduction spin polarization, an^, aii are the partial conductivities of the 
spin-up and spin-down electrons in ith layer, ai — + 0^ is the total 
conductivity of the layer (see [19] for detail). With these assumptions, the 
functions K{t) and P{t) take the form 

P(t) =jHa^^=Po+ Pi COS fit, (6) 
JO 
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P(t) 

= — -T7 =Ko + Ki cos nt, (7) 

OtT^M 

where a is the dimensionless constant of the sd exchange interaction, 

eHaL . . 

Jo = T^, (8) 

fiB is the Bohr magneton, jo has the meaning of the threshold direct 
current density, which corresponds to occurring instability of the antipar- 
allel configuration of the magnetic junction under dominating injection 
mechanism without external magnetic field. 

Let us make the following substitution in Eqs. ([2]), ((3|: 

M^{t) = X{t) exp I ^ sin Q.t 

My(t) ^ Y{t) exp I - ^'^^^^ sin Qt^ . (9) 
The following equations are obtained for functions X{t), Y{t): 

f + + [7(ft + i/.^.)+Po]y 

+KoMzX = -{Pi - KKi)Ycosnt, (10) 

^ - ^ - [7 {h. + HaM. + A-kMTi,^ + Po] X 
+KoM^Y = +{Pi - K,Ki)Xcosnt, (11) 

The exponent in Eq. ([9]) varies within finite limits from exp 
to exp so that stability conditions of the zero solutions for 

X{t), Y{t) coincide with those for Mx, My in the linear approximation 
used here. 

It is seen from Eqs. pop . (jlip that the torque mechanism contribu- 
tion to the effects related with the alternative component of the current 
(Ki) appears with the small factor k and is absent when the damping is 
neglected. 

Let us make the Fourier transformation in Eqs. pO|) . pi|) . The equa- 
tions for the Fourier components with cj frequency take the form 

-iujX{uj) + [7 (^H, + HaM^j + Po - JKcuj y(Lj) 

+Ko'^:,X{ui) = -e [Y{lu + + Y{uj - Q)] , (12) 

~iujY{uj) + (^H:, + HaMz + 47rMM2^ -I- Po - X{uj) 
+KoM.Y{u)) = +e [X{uj + 0.)+ X{u - Q.)] , (13) 

where 

e^'^{P,-^K,). (14) 
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There are "shifted" Fourier components with u ± $7 frequencies in the 
right-hand side of Eqs. (|12p . p3p . If one writes the equations similar 
to pup . (Ilip for such components, then components with oj and u) ± 2Q, 
frequencies appear in the right-hand side of the equations. Continuing 
this process leads to an infinite set of equations coupled via e parameter of 
frequency dimension. We assume this frequency to be low compared to the 
other characteristic frequencies of the problem, namely, the modulation 
frequency Q, and the system natural frequency defined below. Therefore, 
only the component with oj frequency may be kept in the right-hand side 
of the equations for the w ± f2 components, so that closed set of equations 
is obtained. 

In absence of the alternative component of the current (e = 0), the 
frequency w is determined by the zero determinant condition for the ho- 
mogeneous set of equations 



A(aj) = + -2j/3a; + cJo =0, (15) 



where 



2 



+ PoM.) {Q.y + PoM.) + kI, (16) 

(17) 



2 \ K 

= 7 [hJI, +Ha + 4tvM) , Qy=-f (^hJ^, + Ha) 

(usually ttx ^ fiy). 

The fluctuations with frequency determined by Eq. (|15p become un- 
stable even if one of the following conditions fulfils: 

tJo < 0, /3 < 0. (18) 

In presence of the torque mechanism only (Pq = 0), only the second 

of conditions (|18p can be fulfilled at Kq > -{n:, + Qy), = -1 (the 

antiparallel configuration) . In presence of the injection mechanism only 

{Ko = 0), the first condition fulfils at Po > Qy, Mz = — 1. By choosing 
the external magnetic field close enough (in magnitude) to the anisotropy 
field (0 < Ha — Hz "C Ha), the injection mechanism instability threshold 
can be lowered significantly. Such a field does not influence practically the 
torque mechanism threshold because of Q,x S> fiy. Under such conditions 
the torque mechanism may be neglecting, so that Ko = may be laid in 

Eqs. (nni, (HH). 

We seek solution of the dispersion equation in presence of the high- 
frequency component in the form uj = ujq ^ i/, \u\ <C ji^ol- We are in- 
teresting below in stability of the solution under fulfilling the parametric 
resonance condition = 2a;o. 

By virtue of small e parameter, we may neglect the terms in the dis- 
persion equation where that parameter is multiplied by small quantities 
K and V. 

Under actual conditions, AnM ^ Ha, H, Pq inequalities fulfill usually. 
With these inequalities taking into account, the following parametric res- 
onance equation is obtained: 

A(i.j)A(a; - n) = {A-K-yMef. (19) 
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We obtain from Eq. (fT9t 



w = LJo - 27riK7A/ ± (20) 

1^0 



An-yM [7 (^i/.M. + i/.j + PoA/.J . (21) 

The following condition of the parametric instability is obtained from 
Eq. (O: 

e > KCJo, cJo > 0. (22) 
It follows from Eqs. psp and (|22|l that two instability mechanisms, 
namely, a nonparametric one (which is due to the spin polarized direct 
current in absence of modulation) and parametric one due to modulation, 
compete each other, so that only one of them is realized for given config- 
uration (parallel or antiparallel) . However, a situation is possible when 
instability of the antiparallel configuration is created by the nonparamet- 
ric mechanism, while the instability of the parallel configuration is created 
by the parametric one. It is possible, also, that both the coUinear con- 
figurations become unstable due to the parametric instability in absence 
of the nonparametric one. It is necessary to compare the corresponding 
instability thresholds to know what of the variants realizes. 

The frequency ujq takes different values for the parallel (Mz = +1) 
and antiparallel (M2 = — 1) relative orientation of the pinned and free 
layers, these values are denoted below as i^g*'. At \Hz\ < Ha we have 
(4+^) > 0. It means stability of the parallel configuration in absence of 
the current modulation and possibility of its instability under high enough 
amplitude of the high-frequency component. 

From Eq. (pTj) with Eqs. ((6]), ^ and (fTl)) taking into account, we 
obtain the threshold value of the alternative component amplitude j cor- 
responding to occurring parametric instability 

7(±) _ "iei^LJo^^ 

At typical parameter values a ~ 10*, r ~ 10~^^ s, k ~ 10~^, Qi ~ 
10-\ L ~ lO"'' cm, ~ 10^° s-i we have j^'if ' ~ 10'' ^ lO'^ A/cm^ 
The threshold direct current density 

^M1±H,/Ha) (24) 

may be both larger and smaller than j^^'^; this depends, in the first place, 
on the damping constant k and the external magnetic field. Therefore, 
under the mentioned condition j < j a, situation is possible when the 
parametric instability threshold is exceeded, while the nonparametric in- 
stability threshold still does not reached. 

Note that the threshold density of the high-frequency component j^j^^ 
for the antiparallel configuration may be both smaller and larger than 
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Figure 2: Phase diagram for the antiparahel configuration (see explanation in 
the text). 



similar quantity jy^ for the parallel configuration; the latter situation 
takes place at Hz < 0, j/jo < \Hz\/Ha. 

Let us consider the case j^''' < At j < ][^\ j^^'^ < j < Jth' 

the parametric instability of the antiparallel configuration takes place with 
stability of the parallel one; in such a case the initial antiparallel configura- 
tion switches to parallel one. At j < \ j > jih^ parametric instability 
of both coUinear configurations occurs, so that an oscillation regime takes 
place. In the opposite case, j^j^' > the parametric instability of the 
parallel configuration takes place at j < < j < Jtt' and both 

collinear configurations at j < \ j > jth^ ■ 

At j > j'jj ' the nonparametric instability of the antiparallel configu- 
ration occurs (because (^ljq < 0); in this case, as mentioned above, 
parametric instability of this configuration is suppressed. However, at 
j > jl^'' the parametric instability of the parallel configuration takes 
place. 

In Fig. [2] a phase diagram in ^j, coordinates (with dimensionless 

variables Jdc = j/jo, Jac = j/jo) is shown for the antiparallel configura- 
tion at H = and various values of the parameter A = Ak{-kM / HaY^'^ ■ 
The straight line bounds the range of the high-frequency current density 
values according with the condition j < j. On the right of the dashed line 
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Jdc = (th) = Jth /jo the range is placed of the nonparametric instabil- 
ity of the antiparallel configuration where < 0. The parametric 

instability of this configuration at given value of the A parameter takes 
place in the range between that line and the curve corresponding to that 
value of A. 

Thus, instability of any of two collinear configurations, as well as both 
the configurations can be created under amplitude modulation of the spin- 
polarized current. This opens possibility of the switching any of two 
configurations to another (with the same direction of the total current), 
as well creating an oscillatory regime when both collinear configurations 
arc unstable. 

The experimentally observed lowering the threshold of switching by 
direct current in presence of a high-frequency component may be related 
with appearance of the parametric instability creating by that component 
under direct current density lower then the threshold mentioned. 
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